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Abstract. We examine the problem of gauge dependence of the 2PI effective action and its 
^-derivable approximations in gauge theories. The dependence on the gauge- fixing condition is 
obtained. The result shows that 4>-derivable approximations, defined as truncations of the 2PI 
effective action at a certain order, have a controlled gauge dependence, i.e. the gauge dependent 
terms appear at higher order than the truncation order. Furthermore, using the stationary point 
obtained for the approximation to evaluate the complete 2PI effective action boosts the order at 
which the gauge dependent terms appear to twice the order of truncation. We also comment on the 
significance of this controlled gauge dependence. 



I. INTRODUCTION 



Perturbative approaches to the study of equihbrium and non-equilibrium properties of hot and dense media may 
\ lead to inconsistencies and are often plagued with infrared divergences. These problems are linked to the fact that 
calculations in terms of the bare quantities of the underlying quantum field theory (and perturbative approximations 
thereof) fail to describe the collective phenomena in the medium. A strategy to tackle this handicap of the theory 
is to work with dressed quantities, in which the most relevant effects of the interacting ensemble are accounted for. 
These dressed quantities are obtained by means of non-perturbative resummation schemes, which usually involve 
solving a set of self-consistent equations. 

\ An arbitrary resummation scheme will however not guarantee that the conservation laws of the original theory 
are preserved by the dressed quantities. A way to solve this problem is by formulating the scheme in terms of an 
action functional that respects the symmetries of the original theory. A particular kind of such action functionals 
was first introduced in the study of non-relativistic Fermi systems by Luttinger and Ward De Dominicis and 
Martin and Baym and later generalized to relativistic field theories by Cornwall, Jackiw and Toumboulis [Q. 
These functionals, which are derived from the so-called 2PI effective action, involve a diagrammatic expansion in 
terms of two-particle irreducible (2PI) skeleton graphs. A particular choice for an action functional is obtained by 
truncating this diagrammatic series. This defines what is called a ^-derivable approximation. A variational principle 
applied to the resulting action leads to a set of self-consistent equations from which the dressed quantities are obtained. 

A manifest advantage of such a functional formulation is that global symmetries of the original theory are preserved. 
Additionally, the variational principle used to determine the dressed quantities guarantees thermodynamic consistency 
All these useful properties make ^-derivable approximations a very attractive mathematical framework for the 
study of properties of high-energy plasmas. In particular, they may prove useful for QCD plasmas, whose interest 
has grown in recent years due to the possibility of creating quark-gluon plasma in heavy-ion collision experiments at 
Brookhaven and CERN. Calculations of thermodynamical quantities such as the entropy ||] and free energy ||] have 
been achieved using these methods. In those calculations a resummation of the physics encoded in the hard thermal 
loops (HTL) was performed. Important to mention is the fact that, due to the remarkable symmetry properties of the 
HTL, the results in Refs. ^ are manifestly gauge invariant. Non-equilibrium properties can also be formally studied 
within these approximation schemes |^] . They could be used to shed some light on important issues such as thermaliza- 
tion and loss of initial correlations. Very interesting results in this direction have been obtained |^ with scalar models. 

However, an extension of these approximation schemes beyond the HTL regime in the study of QCD plasmas is 
still lacking. There are two main problems involved. One is that renormalization seems to be a non-trivial issue, as 
shown in explicit calculations for scalar theories ||^. To deal with this obstacle, a recent approach based on BPHZ 
renormalization has been proposed by van Hees and Knoll 0. The other main problem is the fact that gauge 
invariance may be lost in the approximations. This is because, in general, the solutions for the dressed propagators 
and/or vertices do not satisfy Ward identities. In particular this implies that thermodynamical quantities computed 
within these approximations will suffer from gauge dependence. This pathology shows up as an explicit dependence 
on the choice of gauge-condition. 
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In this paper we study the problem of gauge dependence of the 2PI effective action and its ^-derivable approxima- 
tions. In Sec. II we review the general formalism of <i>-derivable approximations and introduce the notation to be used. 
In Sec. Ill we apply the formalism to gauge theories and determine the dependence of the 2PI effective action under 
a change of the gauge-fixing condition. From the result one sees that the 2PI effective action is gauge independent at 
its stationary point. This was already shown for the IPI effective action and expected from general arguments pT[ . 
In Sec. IV we apply the result of Sec. Ill to the ^-derivable approximations that result from truncating the 2PI 
effective action at a certain order. We show that these approximations have a controlled gauge-fixing dependence, i.e. 
the gauge-dependent terms appear at higher order. We discuss in Sec. ^ that the use of <I>-derivable approximations 
restricts the choices of gauge fixing available, if they are indeed to be good approximations to the exact theory. This 
prevents the high-order gauge-dependent terms to take arbitrarily large values, such that the gauge dependence will 
be indeed controlled. 



II. 2PI EFFECTIVE ACTION AND ^-DERIVABLE APPROXIMATIONS 



The generating functional for correlation functions can be written as 

Z[J,K] = j 2?^e*{^['^l+'^-'^'+2'^'-f^'^'^'}, (1) 

where S[(f\ is the action, represents the fields and the J and K are auxiliary external sources. We use a shorthand 
notation where Latin indices stand for all field and current attributes (i.e. f{x) — > ipi) and summation and/or 
integration over repeated indices is understood, i.e. JiLp^ = j d'^x J{x)f{x)} . The generating functional of connected 
diagrams W is defined from Z as 



W[J,K] = ^i\og{Z[J,K]). 
The expectation value of a functional 0[lp\ is given by 

(C 03 ) = p ■ ^ = lO 



5{iJ) 



W. 



(2) 



(3) 



Mean fields 0* and connected correlation functions G"^^ can then be obtained by functional differentiations of W[J] 
as 



. 5W 



S{iJi)S{iJj 



= G'\ i 



5{iJi)5{iJj)5{iJk) . 



(4) 



Functional differentiations of W[J, K] with respect to the bilocal currents K may generate also disconnected dia- 
grams. For example, differentiating once with respect to K leads to 



. 6W[J,K] 



1 



(5) 



A functional Legendre transform in the mean field 0* and the two-point function G*-' leads to the so-called 2PI effective 
action 



,G] = W\J,K] - J,, 



2 ^ 



iV^' + G'^) , 



From its definition one can derive the relations 

smG] 



^~J,~K,j(ty> and ^^^^TT^ = "^^y- 



(6) 



(7) 



The time integration involved in this functional product can also run along a contour C in the complex plane such as the ones used in 
the real and imaginary time formalisms of thermal field theory. This detail will however not be important in our calculations, so we will 
omit the subscript C in the integrations. 
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With the help of Eq. (0) one can write the expression for the expectation value of a functional 0[ip] in terms of the 
2PI effective action as 

{0[ip]) = e-*r[0,G] J 

The 2PI effective action can be cast into the very convenient form ^ ^, ^, ^, |^ 

r[0, G] = ^o[0] + «cTr {log (G-i) + G (Gp ^ - G-^)} - i<i>[(/., G]. (9) 

where 5*0 is the free part of action, Gq is the bare two-point function (—iS'^ So[lp]/ Stpd(f) ^ and c is a constant equal 
to 1/2 for bosons and —1 for fermions. The functional G] consists on the sum of all two-particle-irreducible (2PI) 
skeleton diagrams with bare vertices and dressed propagators. In this context skeleton diagrams are those without 
self-energy insertions. Non-2PI diagrams with mean field insertions are also included in the definition of <i>.^ For 
example, in the case of a theory with quartic interactions (such as A^''), and using Dyson relation G^^ = Gq^ + ill 
between the two-point function and the self-energy 11 , the above expression can be written graphically as 

r|*,G1^5„|*| + ic{Q.(^)-,{±;^.|^.|g.J;a0..±^.., (10) 

where the thick lines are dressed propagators, the small lollipops are the mean fields and the cross-hatched blob is 
the one-particle irreducible self-energy 11. Writing the effective action in terms of skeleton diagrams makes possible 
to incorporate higher-order effects into the propagators with correct prefactors. 

In practice one has to restrict to an approximated version of the 2PI effective action that results from considering only 
a certain subset of diagrams in the functional <&. This defines a ^-derivable approximation. One typically considers 
the loop expansion of skeleton diagrams as pictured in Eq. ( p^ ) (and refered to as the skeleton-loop expansion in 
the following), and truncates it at a given order. In this way the 2PI effective action F is split into two pieces: the 
truncated part Fq, and the higher order part Fi. Then one takes Fq as the approximated effective action. This action 
defines approximate mean fields and two-point functions 0ap and Gap, which result from the stationarity condition, 
i.e. from the implicit functional equation Wj for vanishing sources J and K, as 



sro[(P,G] 



0ap : Gap 



and 



<5Fo[0,G] 



SG 



0ap ! Gap 



0. 



(11) 



These approximate mean fields and two-point functions defined from the truncated action Fg differ from the exact 
ones, which are obtained from the stationary point of the complete 2PI effective action as 



5r[0,G] _ <5(Fo-f Fi) 



and 



5T[^,G] <5(Fo + Fi) 



6G 



SG 



0. 



(12) 



We end this section by noting that one could also construct more general effective actions by including higher-point 
external sources Ji,Kij,Lijk ■ ■ ■ into the functional W = W[Ji, Kij, Lijk, ■ ■ ■] and performing a Legendre transform as 
follows 



T[(j)i,Gij,Gijk, ■ ■ ■] — W[Ji,Kij,Lijk, 

1 



-J,: I 



2^- 



V^' + G'') - -Ujk (G'^'= + 3G'>^ + (/'W'^) - . . . . (13) 



This form of the effective action can be rewritten as a diagrammatic series in terms of skeleton diagrams of the n-point 
vertex functions and can be used for generalized ^-derivable approximations. 



^ In the literature <& is usually defined in such a way that it only contains strict 2PI diagrams. This involves a redefinition of the 
action to include mean fields and tadpoles. We prefer the above notation where all interaction parts are placed in Of course, both 
definitions agree when (/> = 0. 
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III. GAUGE-FIXING DEPENDENCE OF THE 2PI EFFECTIVE ACTION 



We consider now the case of a pure Yang-Mills theory with gauge group G. Its action is given by 

Sym = ^J d'x ^F^,{x)Fr{x), (14) 

where i^^^ = F^^ Ta ~ d^A^ — d^A^^ — g [A^, Ay] is the field-strength tensor of the gauge field A^ — A^^^Ta, g is the 
(unrenormalized) coupling constant and Ta are the generators of the Lie algebra of the gauge group G. 



The action is invariant under gauge transformations U{x) G G of the gauge potential A^ 



A^ ^ = U{x)A,,{x)U-'{x) - [d^U{x)] U-\x). (15) 

This invariance implies that the functional integrals over gauge field configurations are ill-defined. One gets around 
this difficulty by the Faddeev-Popov gauge-fixing procedure, which introduces a gauge-breaking term S'gf into the 
action. In the context of BRS-quantization this term is realized in a useful manner by introducing some auxiliary 
fields: the Faddeev-Popov fermionic ghost fields c° and and the bosonic Lautrup-Nakanishi fields B". The gauge- 
fixing is implemented through the condition C°[yl] = 0, where a typical choice is the covariant gauge C"'[^] — 
The gauge-fixed action then reads 

S = Sym + Sgf = J d^x ^-\F%{x)Fr{x) - Ux)^-^^{DA^)), + Ba{x)C'^[A] - i^i3„(x)i?"(x)| , (16) 
where Z)^ = dp,, — igT'^Aa is the covariant derivative and ^ is the gauge-fixing parameter. 



The action obtained by adding this gauge- fixing term is no longer invariant under local gauge transformations (15). 
However, it is invariant under BRS transformations, which are defined as 

S^^sA; = e{DA\ 
f^BRsc" = if-gc^, 

S^usB'^ = 0, (17) 

where e is an infinitesimal global anti-commuting parameter and is a short-hand notation for {Tac"'){Ti,c^) ~ 
l/2[r°, r'']c"c''. The Lautrup-Nakanishi field B has been introduced to ensure the nilpotency of the BRS charge 
Qbrs, defined as (5brs — cQbrs- It allows for a convenient rewriting of the gauge-breaking term as a complete BRS 
variation 

.4. n 



^GF = QbusJ d^x |-Cc„(x)B'^(x) - c,(x)C°[A]| = geos*. (18) 
Using the notation of Sec. II the generating functional Z[J, K] for this gauge theory can be compactly written as 

Z[J,K] [ X'(pe'{'^™+'3™s*+'^iv'-H5¥''^i3v'} = gj{r[0,G]-HJi0'+iKij(0>^+G")}^ ^g-j 



where (f denotes collectively all fields {A'^{x),c°'{x),c°'{x),B°'{x)}, J and K denote all their associated currents 
{Ja, Jc-, Jc, Jbt KaAt Kcc, Kbb}, and Latin indices stand for both space-time and group indices, i.e. A'^{x) —> A^. 
This notation is used to allow one to write formulas in a compact way, though one should bear in mind that the 
ghost fields c and c and their associated local currents Jc and Jg are anti-commuting variables. is a ^-dependent 
infinite constant generated during the Faddeev-Popov gauge-fixing procedure. Its gauge parameter dependence can 
be seen already in the free theory and can be absorbed into the action by rescaling the ghost fields by Hence 
this constant will not play a role in the following. 

Having set up the notation we turn now to study the gauge dependence of the 2PI effective action. We study 
how it transforms both under a change of the gauge-fixing condition [A] C° [A] + AC" [A] and gauge parameter 
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^ ^ ^ + A^, or more generally, under a change ^ + A^". Under this shift of gauge condition, the effective action, 
the mean field and the two-point function respectively change as 

r ^ r' = r + Ar , -> 0' = + A0 and g ^ g" = g + ag. (20) 

The currents J,; and Kij are taken to be gauge independent since they are external. This fact allows us to calculate 
immediately from Eq. (|7[) how much the first functional derivatives of the effective action vary under the gauge-fixing 
change, obtaining 

, Mr \ . fsr\ ^ ST 



The first functional derivatives of F are used to find the stationary point by setting them to zero as done in Eq. (12). 
From Eq. (|2^) one notices that the stationarity condition itself is only gauge invariant if it is realized simultaneously 
for both arguments (f) and G. 

To compute the variation of the effective action F itself one can use the relations (|^) to cast Eq. ( |l9| ) into the 
convenient form 

^r^l4>.,G] ^ J -p^^'^{SYM + QBRS'Sf-{V-<t>)^-i^-l{V^-'t>^){Vj-<pJ)-GiJ]J§--] ^22) 

For simplicity and later convenience we denote the field combinations (93 — and [{ip — 4')ii'P ^ ^ Gij] that 
appear in the exponent of Eq. (|2^) as ipi and G^ respectively. These have the property that their expectation values 
(ifii) and (Gij) vanish. 



After a change of gauge condition ^' ^ ^P' = ^' + A^" equation (^2|) becomes 

g^r' ^ g»(r+Ar) ^ y |g'{^YM+QBns*-¥'.i^-G,,7%} ^^IqbesA^+a^.^+ac.^+a^.a^,^} |^ ^23) 

which using the notation of Eq. (||) leads to 

giAr ^ ^^i{QBRsA*+A0.^+AG.,^+A0,A0,^}^ ^^4) 

This result is valid for any finite change in the gauge- fixing conditions. To proceed further we restrict ourselves to 
infinitesimal variations A^*. Then one can expand both sides of Eq. (24) to obtain 

AF = (gens A*) + A0,;|^^ + AGy + ©(A^) (25) 

where we used the fact that Acj) and AG are of order 0{A'$). This can be easily checked. Indeed, following the same 
steps as to obtain Eq. (24) one gets for the mean field 

which using Eq. (U) reduces to 

^ + ^^- (g.gBBsA.)^ (27) 

and expanding in A^* yields 

A0, =i (^.QbbsA*) +0(A2). (28) 
Similarly one computes the variation of the two-point function obtaining 

AGy = ^ (GyQBR.sA*\ + ©(A^), (29) 
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which verify our statement above. Moreover, Eqs. ( |28D and (|29|) can be used to write Eq. (^5[) as 

Ar = (QbbsA^') + z ((p.QbrsA*) — + i (G,,QBnsA* ) -— + 0(A2). (30) 

One expects the stationary point of the effective action^ i.e. when its functional derivatives are set to zero, to be 
gauge-independent. That is still not obvious from Eq. (pQ) since it appears that the first term in the r.h.s. would 
not vanish. For that, one can make use of the following trick [ p^ . Consider the expectation value of the gauge-fixing 
change A 4', namely 

(A*) = e-'^ j 2?<^A^eH''™+'^™^*-^'^-^-^}. (31) 

One can do a BRS transformation ip Lp + eQoBs'/' on the field variables in the path-integral. This transformation 
leaves the measure invariant and only amounts to a shift of the integration variable, so the equation remains the same. 
The l.h.s can be however rewritten so that 

(A*) -e-''!' j Vp (A* + eQ^„^A*)eH^™+'2™'*"^"^"^'^'^"''^™'^'^"''^™'^'^'^ (32) 

where we have used the fact that lS}i![ip + eQ-a-as^] — A^* -I- eQensA^', which follows from the definition (|l^). The BRS 
charge Qbrs appearing in this expression does not operate on the mean fields (j) and two-point functions G that are 
part of if> and G, but only on the fields to be path-integrated over. Expanding the r.h.s. of ( ^2|) in the anticommuting 
parameter e leads to 

(QensA^') - -I <^A* Qens (^^^) ) " ^ Qbhs (^U ^) ) ' (33) 

where the quantities have been reorganized so that the equation is valid for all fields, both commuting (yl^ and B) 
and anticommuting (c and c). Combinations like ip6T/6(j} or GST/SG are always commuting so it is preferable to 
have them in this form. 

The same procedure can be applied also to the expectation values (^aA^) and (GabA^) to obtain 

(paQsn^^'f) = (A^'QB„s<?a>-j<^^aA^'QBns(^^j^^^-«^^aA^'QBns(^Gjfe^^^^, (34) 
Gab QbbsA*) = (a* QBBsGah) - i (CabA-^ Qbbs (^J^) ) ~ * (^afc A* Qbbs (^j'c^) ^ • (35) 

These results enable us to write the ch ang e in the effective action AF only in terms proportional to its functional 
derivatives simply by substituting Eqs. (33)-(^5|) into Eq. (|3^). One notices that the first terms coming from the r.h.s. 
of Eqs. ( ^4| ) and (|3^) cancel exactly those that come from Eq. (33) when both are substituted into Eq. (|30|). In this 
way terms with a single functional derivative do not appear in AF. After some rearrangements one is then left with 

AF ^ 1 (avI^Qb. (^^^^.|-)) + (a*Qbbs JJ;)) + I (a^Qbbs (g., J-G. J-)) +G(A^), 



(36) 



which can be cast into the compact result 



AF[</),G] = i/A*QBn,s(^^.^ + G,fe^) )+0(A^), (37) 



where the average and Qbrs only apply to the fields contained in ip and G. 



Eq. ( |37| ) gives the variation of the 2PI effective action caused by a change in the gauge condition and is the main 
result of this paper. One sees that when the functional derivatives of F are set to zero this variation vanishes, and 
then the effective action is gauge-fixing independent. This situation occurs precisely at the stationary point, i.e. at 
the exact mean fields 0ox and two-point functions Gcx- 



7 



However, the quantities 0cx and Gox are not gauge-fixing independent themselves. Indeed, one can explicitly 
compute their gauge dependence by applying the condition ([l^ ) to Eqs. (28) and (p9|), obtaining in this manner 

AcPl^ = * (Avl/g,„s (^cx)') , (38) 

AG^i = i<^A*Qbrs(Gcx)''^. (39) 

For the case of the effective action F [(jji , , G^-fe , ■ • ■ ] including higher-point correlation functions, the same procedure 
leads to the generalized result 

AF ^ -i (^A^ Q.. (^.^ + G., J- + G.....^ + . . y ^ + G(A^). (40) 

where the quantities Gijk... are given by 

Gi = = - (41a) 
Gy = ((^-0),((^-0), -G,,, (41b) 

GijTc = ((^ - 0)4<<2 - 0)j('^ -</»)*;- Gijfc - G.ij(<^ - 0)a: - GjA;((;£' - (j!))i - GA:i((^ - (?!))j, (41c) 

Gijki = ("^ - 0)i(<<5 - (t>)j{v - 4')k{v - 0)i - Gijkiv - 't>)i - Gjki{(p - 0)1 - Gfcii((y9 - - Guj{(p - 0)^ - 

GijGki - GifeGj/ - GijGfcj - {G^j [{(p - (j))k{^ - (j))i - Gui] + 5 permutations } - Gijku (41d) 

etcetera. 



IV. GAUGE-FIXING DEPENDENCE OF ^-DERIVABLE APPROXIMATIONS 

Our main interest is to study the gauge dependence of <i>-derivable approximations. As previously mentioned, 
they are obtained after one truncates the skeleton-loop expansion of the 2PI effective action at a certain order. For 
definiteness let us consider a truncation at L loops, which translates into a truncation at 0{g^^~^) for the coupling 
constant g. Then F is split into two pieces 

F[0, G] = Fo[0, G] (of 0(g2^-2)) ^ 0(5'^)), (42) 

where the truncated part Fq is used to generate approximate mean fields (/)ap and two-point functions Gap from the 
stationarity condition (|l|). 



This splitting of the effective action can be performed directly on the result (^7|) for the variation AF under a shift 
of gauge evaluated at the approximate mean fields (/)ap and two-point functions Gap 



1 / A T ^ /(5Fl 



5Fi 



A(Fo+Fi)[(/.ap,Gap] = --(A*Qens(xr . ^ ^^.-P + IT^. ^ G,fe, ap ) )+0(A^), (43) 



Gap ' SGjk 



</>ap,Gi 



ap 



2\ 



where we used the fact that 0ap and Gap correspond to the stationary point of Fq 



On one hand, Eq. ( [4 3D implies that the truncated effective action Fq evaluated at its corresponding physical mean 
fields (/)ap and propagators Gap is gauge independent up to the order of truncation, i.e. 0{g^^~^). This is so since Fi 



is of order 0{g^^) and the r.h.s. of Eq. (43) is of order O(F^), so to first order AFq w -AFi w Gig ) 



On the other hand, Eq. ( [13[ ) tells us that the complete action F evaluated at the approximate mean fields and 
propagators obtained from Fq is gauge-fixing independent up to order 0{g^'"), i.e. twice the order of Fi. This is a 



^ Since one works at the stationary point of Fq instead of the exact one obtained from F this imphes immediately from Eq. (pt that 
expectation values (• • • ) are here evaluated at the values of the currents given by J' = —5Ti/5(f) + 24iSTi/SG and K' = —2 STi/5G. 
However, since Fi ~ 0{g^^), the expectation values are in first approximation equal to those obtained with vanishing currents. 
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consequence of having the square of the functional derivatives of Fi in the r.h.s. of Eq. ( ji^ ) and can be understood 
with a diagrammatic argument. To see that, first note^ that the diagrams in the loop expansion of the 2PI effective 
action F are skeleton diagrams, so without any self-energy insertions. The approximate propagator Gap obtained 
from truncating the skeleton series of F to Fq at L loops (or at 0{g^^~^)) is the solution of the variational condition 
([T^), which can be interpreted as a dressing of the bare propagator with all the self-energy contributions that come 
from cutting one line in the 2PI diagrams of Fq. Evaluating the effective action F at Gap entails substituting this 
propagator in the diagrams of the skeleton-loop expansion. The outcome can be expanded perturbatively to compare 
directly with the usual perturbative loop expansion of the IPI effective action. One can check that both expansions 
match perfectly up to 2L loops, or 0{g'^^~'^). They differ at 2L -I- 1 loops because, by construction, diagrams that 
would result from dressing skeleton diagrams of Fi with self-energy contributions to Gap coming also of Fi, do 
not appear in the expansion of the skeleton series considered. However, they are present in the perturbative loop 
expansion. The importance of the fact that both expansions match up to 2L loops is that, since the perturbative loop 
expansion is gauge-invariant at every loop order, one can immediately conclude that so must be the skeleton-loop 
expansion of F[Gap] up to 2L loops, or, in other words, up to 0{g'^^). 

In this manner, Eq. (^3|) shows that <i>-derivable approximations, as truncations to the 2PI effective action, have a 
controlled gauge- fixing dependence, in the sense that gauge dependent terms appear at higher orders. 



V. CHOICE OF GAUGE CONDITION 



A large body of experience with gauge theory has led to the common view that one should not tamper with gauge 
invariance. Yet, we explore here the possibility of accepting a controlled amount of gauge dependence in the compu- 
tation of physical quantities. The question is then, what is a good choice of gauge fixing? To be specific, consider the 
class of covariant gauges described by G° = d^AI^^. Then we have to decide on a reasonable choice for the gauge pa- 
rameter ^. Evidently, ^ should be such that it does not upset the assumption that Fi may be neglected compared to Fq. 

That such upset can happen is easier to see in the more familiar perturbative case. There we have the loop 
expansion in terms of bare propagators. Consider for simplicity diagrams without ghosts. The gauge propagators 
have a longitudinal part proportional to ^, so in a given diagram with / internal lines we would have the factor 
In terms of the number of three- and four-point vertices V3 and V4 and the number of loops L it can be written 
as ^2L-2+V3/2^ Together with the powers g^^^^ in the bare coupling constant g, the diagram has an overall factor 
■ Taking <^ big enough, say |^| > 1/g, various terms belonging to different orders of 5 in the perturbation 
expansion would be shuffled. This will evidently upset our ordering principle. 

In a <I>-derivable approximation, however, we consider the loop expansion in terms of dressed propagators where 
their ^-dependence is not clear a priori. For that one needs to find the stationary point of Fq. And this is done after 
assuming that Fi may be neglected compared to Fq. Provided we had the explicit form of the dressed propagator, 
an argument similar to the one above would give the range of ^ that is allowed without upsetting this assumption. 
Unfortunately, finding the dressed propagators is in general a formidable task. We nevertheless venture the following 
argument that a good choice for is in the interval (0, 2). 

Assume that the ^-derivable approximation gives indeed an approximation to the path integral 

Z = jvAVcVcexp -J J d^x f^^F^^, + cd'^id,, - iA^)c+ ^{d'^A^f'^ , (44) 

where we have integrated out the B field and rescaled A A/g, (c, c) {c/y/g, c/y/g). Then g^ in the above action 
is also the ordering parameter in the skeleton- loop expansion of <&. If we do not want to upset this power counting, 
should be treated of order one as g^ 0. For finite g^ it seems best to choose 1/2^ of the same magnitude as the 
other numerical coefficients in the action, which are 1/4 and 1/2 for F^, and 1 for the ghost terms. So this suggests 
the choice ^ in the range 1 — 2. Saddle point arguments for 5^ — > are not upset by letting also ^ ^ 0, so it is 
reasonable to allow also values of ^ 0. On the other hand, ^ ^ 1 would upset the longitudinal parts of the saddle 



We take = for simplicity. 
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point regions in functional space (as for the perturbative case above) . To allow a continuation of the path integral to 
imaginary time, ^ has to be positive. All these arguments lead us to conclude that ^ is best taken in the range 0—2. 

VI. CONCLUSIONS AND COMMENTS 

In this paper, the gauge dependence of the 2PI effective action that defines a <i>-derivable approximation of a 
gauge theory has been determined. To obtain it we used its definition as a Legendre transform of a generating 
functional with bilocal sources and the BRS symmetry of the underlying Yang-Mills action. As expected on general 
grounds, the result shows that the 2PI effective action is gauge independent at its stationary point. Furthermore, the 
result has been applied to study the gauge dependence of <i>-derivable approximations, defined as truncations of the 
2PI effective action at a certain loop order. Even though correlation functions derived within these approximation 
schemes are known not to fulfill the Ward identities required by the gauge symmetry, it has been shown that the 
truncated effective action defined at its stationary point has a controlled gauge-fixing dependence, i.e. the explicit 
gauge dependent terms appear at higher order. Furthermore, if one uses the stationary quantities of the truncated ac- 
tion to evaluate the complete 2PI effective action, the gauge dependence appears then at twice the order of truncation. 

These features might be interesting for the computation of thermodynamical quantities derived from the 2PI 
effective action in gauge theories, such as the pressure and entropy. The authors of |^ have calculated the entropy of 
the quark-gluon plasma with an approximate ^-derivable approximation, but in order to achieve gauge independence, 
they had to sacrifice the self-consistency guaranteed by working at the stationary point, hence the word approximate. 
Their approach was nevertheless strongly motivated from a quasiparticle picture of the quark-gluon plasma, which can 
be used Q to describe the lattice results In any case, our considerations suggest that <i>-derivable approximations 
may allow a systematic method for computing thermodynamic functions without having to sacrifice its remarkable 
properties. Gauge-fixing dependent artifacts would appear at high orders, thus making the approximation controllable. 

It might seem unsatisfactory that the gauge dependence is not completely removed in ^-derivable approximations. 
Yet, we propose here to accept a controlled amount of gauge dependence in physical quantities. We argued that 
<I>-derivable approximations, provided they are indeed an approximation to the exact gauge theory, implicitly restrict 
the choices of gauges available. This prevents the high-order gauge artifacts to take arbitrary values that could 
render any computed quantity physically meaningless. Both the fact that gauge dependent terms appear at higher 
orders and that they are constrained by this restriction makes the error introduced by breaking gauge invariance 
controllable and indicate that ^-derivable approximations may indeed give reasonable answers to physical quantities. 
A detailed examination to quantify those gauge dependent terms would involve solving a ^-derivable approximation 
for gauge theories. So far, the complete solutions for QED and pure glue QCD even at lowest order (2-loop) have not 
appeared in the literature. 

We would like to note that in the derivation presented in this paper we did not discuss aspects related to regular- 
ization and renormalization. This makes the calculations heuristic in some respects. We compared the path integral 
with a skeleton expansion of the effective action, neither of which has been clearly defined. To make the path integral 
well defined a regularization is needed, preferably non-pcrturbative, and this should be compatible with the BRS 
invariance used. With a lattice regularization this is non-trivial [l6| . Another point is that the regularization depen- 
dence needs to be removed, or at least shown to be negligible (in the case of QED 'triviality' is expected to occur). 
Renormalization is a non-trivial issue in <I>-derivable approximations A general renormalization procedure, such as 
the one recently proposed by van Hees and Knoll |jT^ , would be needed in order to have a well defined path- integral. 
A detailed study of these issues in <I>-derivable approximations of gauge theories constitutes the subject of further 
investigations. 
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